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Abstract. The purpose of this article is to develop a theory behind the occurrence of 
"path-integral" kernels in the study of extended determinantal point processes and non- 
intersecting line ensembles. Our first result shows how determinants involving such kernels 
arise naturally in studying ratios of partition functions and expectations of multiplica- 
tive functionals for ensembles of non-intersecting paths on weighted graphs. Our second 
result shows how Fredholm determinants with extended kernels (as arise in the study of 
extended determinantal point processes such as the Airy2 process) are equal to Fredholm 
determinants with path-integral kernels. We also show how the second result applies to a 
number of examples including the stationary (GUE) Dyson Brownian motion, the Airy2 
process, the Pearcey process, the Airyi and Airy2->i processes, and Markov processes on 
partitions related to the z-measures. 

1. Introduction 

The Airy2 process is a universal scaling limit of a wide variety of probabilistic systems 
including random matrix theory, random growth processes, interacting particle systems and 
directed polymers in random media (see [24,41] and references therein). Denoted Airy 2 (-), 
it is defined via its consistent finite dimensional distributions: for t\ < Ct2 <C * * * <C t n ^ 

pfQ {Aiiy 2 (*0 < = det(J - X^2 Xt )L2 {{tl ,..., tn}xM , At ). (1) 

Here x is an operator which acts on functions /: {t\, . . . ,t n } x R — > R as 

Xf{U,x) := l x > s J(ti,x). 

The operator Iff acts as 

n 

Kl xt f(t h x) :=V / dyK^fax^Mfay) 

where K2 Xt (s,x;t,y) is the "extended" Airy 2 kernel given by 

_ (j oo d\e-^Ai(x + \)Ai(y + \) if s > t, 
2 1 ' ' hyh ~ \-/_° oo dAe- A (- t )Ai(x + A)Ai(y + A) if a < t, 

with Ai(x) the classical Airy function. The right-hand side of (1) is the Fredholm determi- 
nant of the identity minus a trace class operator (see Section 3.1 for definition and details) 
and the measure fj, appearing there is the product of counting measure on {t±, . . . , t n } and 
Lebesgue measure on R. 

The formula given in (1) for the finite dimensional distributions of the Airy 2 process 
becomes increasingly cumbersome as n increases. This is due to the n-dependence in the 
L 2 space on which the operators act. When taking a limit of a sequence of operators, or 
their determinants, it is convenient to have the operators all act on the same L 2 space, 
rather than a sequence of different spaces. 
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In Prahofer and Spohn's initial work on the Airy2 process (see Section 5 of [36] for n = 2 
or [22, 37, 40] for n > 2) the extended kernel formula is shown to be equivalent to the 
following "path-integral" kernel formula: 

pf P {Aiiy 2 (ti) < Si }j = det(I-K 2 +P sl e^-^ H P S2 ■ ■ ■ e^-^ H P Sn e^-^ H K 2 ) L 2 {R) . 

(2) 

Here K 2 (x, y) = K^iO, x 'i 0> u) ls called the Airy2 kernel, P s g{x) = l x < s g(x) is a projection 
operator and H = —A + x is called the Airy Hamiltonian because HAi(- — s) = sAi(- — s) 
(A is the Laplacian on R). The dependence on n has been absorbed into the operator 
rather than the I? space, and it is now plausible to take a large n limit. 

The reason we call this a path-integral kernel is because a portion of it can be written 
in terms of the expectation of a certain path-integral. By the Feynman-Kac formula (cf. 
[28]), 



P sie (ti-t2)Hp s2 . . . e {t n -^t n )Hp sJ . {x) = E 



lb 

f(b(t n ))e-^ b{s)ds Hl 



b(ti)<Si 



i=l 



(3) 



where b: [ti,t n ] — > R is the trajectory of a Brownian motion with diffusion coefficient 2 
starting at b(ti) = x. 

Let t\ < ■■■ < t n fill out the interval [i,r] and let Sj = h(ti) for some function h : 
[£, r] — > R. Then as n goes to infinity, the above operator has a limit in trace norm (see 
[22] or Section 4.2 below for details) given by 

r£ r /(s) = E b[e)=x [f(b(r))e- XT b ^ ds l b < h 

where {b < h} denotes the event {b(s) < h(s), Vs G [^ 5 r ]}- Thus, it is shown in [22] 
Theorem 2 or equation (41) below that 

P(Airy 2 (s) < h(s), Vs 6 [€, r]) = det(7 - + r^ e ( r -^ H K 2 ) L2(R) . (4) 

The above formula proved useful in [22] in providing a direct proof that the value of the 
maximum of the Airy 2 process minus a parabola is distributed according to the (GOE) 
Tracy- Widom distribution; and in [32] in computing the joint distribution for the value 
and location (in t) of the maximum. The two-time path-integral kernel formula in [36] 
was utilized to compute asymptotics of the two-time covariance of the Airy 2 process, since 
the extended kernel does not easily yield this. Note that the left-hand side in the last 
formula presupposes the existence of a continuous version of the Airy 2 process. This was 
first shown to exist in [26]. 

What is remarkable about formula (4) is that the right-hand side is simple (despite the 
cumbersome finite dimensional distributions given above) and the event in question in the 
left-hand side has a clear translation into the operator T^ r . As a further application of 
the Feynman-Kac formula as well as the Cameron-Martin-Girsanov formula (see [22]), the 
integral kernel of r^ r can be expressed as 

r e, r ( x ^y) = ^b(e)=x-ez,b(r)=y-r 2 { b ( s ) < K s ) - s 2 for a11 s 6 i^ r ]) . 

where b is now a Brownian bridge run from x — I 2 at time t to y — r 2 at time r (this means 
that F^ r f(x) = f dyTg r (x,y)f(y)). In other words, the probability that the Airy 2 process 
hits a function h can be expressed as the Fredholm determinant of an operator which is 
partly expressed by the probability that a Brownian bridge hits the same function (minus 
a parabola). 

We will now see how these formulas for the Airy 2 process are part of a more general 
result. 
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1.1. Extended kernels and the path-integral kernel in a general setting. There 
are many other examples of extended determinantal point processes (some given in Section 
4) and our aim is to find path- integral kernel formulas for these other processes. This may 
have further applications, although we do not address them here. For example, besides the 
previous work of [37], in [40] a path-integral kernel formula was discovered for the Airyi 
process and used to prove existence of a continuous version of the process and its Holder 
regularity. 

When the Airy2 process was introduced, it arose as the top layer of the multi-layer Airy2 
process, which will be denoted {Airy 2 (i;i) : i £ Z>i,t G R} and is such that Airy 2 (z;i) > 
Airy 2 (j;i) for i < j. As t varies, YlZi ^Airy 2 (i;t) forms an R-indexed collection of point 
processes which has the structure of an (extended) determinantal point process (see [7] and 
references therein) with correlation kernel K^ Lt {s,x]t,y). 

This has the following consequence. For t\ < ■ ■ ■ < t n fix : R — > R and let q t . = 1 — q t . . 
For g: {ti, . . . , t n } — > R define q(g) := J\i=i Qu{g{U)) and likewise define q(g). Then (given 
some conditions on q to ensure convergence - see Section 4.2 below) 



where Qf(U,x) := qt i (x)f(ti,x). The left-hand side above is referred to here as the ex- 
pectation of a "multiplicative functional" of the multi-layer process. When q k (x) = l x >si 
(and hence qt^x) = l x < Si ), the above formula reduces to (1) with Q = X- 

Our first result shows that the type of identity between extended and path-integral kernel 
Fredholm determinants which one gets by equating the right-hand sides of (1) and (2) is 
quite general and dependent on a few structural properties of the kernels. 

In stating our results presently, we leave out a number of technical assumptions (see 
Section 3 for these details) and assume that we have the following collection of operators 
on functions / : {t±, . . . , t n } — > R: 

• For each 1 < i,j < n, Wti,tj (with the convention Wti,u = I)] 

• For each 1 < i < n, K^; 

• A diagonal operator Q such that Qf(U,-) := QufiU-,-) where for 1 < i < n and 
5 :R^R, Q u g(x) := q u (x)g(x). 

Theorem 1.1 (Theorem 3.3, with technical assumptions suppressed). Assume that for all 
l<"i<.7<^< n the following holds: 

• Right-invertibility: W^^t^t^uK^ = K u ; 

• Semigroup property: W^Wt^ = W ti ,t fc ; 

• Reversibility relation: Wt u t ] K tj =K ti Wt i ,t j - 



det(/-QiT xt ) L2({ti jMxR)M) = det{l-K tl +Q tl W tl ,t 2 Q t2 ■ ■ ■ Wt^ ,t n Qt n Wt n ^K tl ) ^ 



oo 



E J] g(Airy 2 (i; •)) = det(J - ({*!,., t„}xR,/x) 



(5) 



.i=i 



Then 



where 



K cxt (t i: x; tj,y) 



W k , tj K tj (x,y) ifi>j, 
-W ti , tj (I - K tj )(x,y) ifi<j, 



and Q t . = I — Q ti ■ 



This is proved in Section 4.2 essentially via linear algebra. 
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Letting Wt u tj = e~^ t i~ ti ' H , K t . = K 2 and Q Si = P Si we recover the equality (2) for 
the Airy2 process. More generally, the result implies that the expectation of a multiplica- 
tive functional of the multi-layer Airy2 process (5) can be expressed in a similar way, by 
replacing each P Si by Q Si on the right-hand side of (2). 

In Section 4 we apply this theorem to a variety of examples of extended determinantal 
point processes such as the stationary (GUE) Dyson Brownian motion, the Airy2 process, 
the Pearcey process, and Markov processes on partitions related to the z-measures. We 
also show how the identity applies to signed extended determinantal point processes such 
as the Airyi and Airy2_s.i processes. In the case of the stationary (GUE) Dyson Brownian 
motion and the Airy2 process, we also obtain the continuum limits of the corresponding 
path- integral kernel formulas (which is likely doable in other cases as well). 

1.2. Ensembles of non- intersecting paths and the path-integral kernel. The multi- 
layer Airy2 process arises as the scaling limit of a variety of ensembles of non-intersecting 
paths (see for instance [27]). The occurrence of extended kernel determinants in such 
ensembles is a consequence of the Eynard-Mehta theorem, which implies the existence of 
an extended determinantal point process structure [17,23,26,33,45]. The equivalence 
of the extended kernel determinant formula with the path-integral kernel determinant 
formula which is given in Theorem 1.1 (see also Theorem 3.3) is via linear algebra, but 
does not indicate why such a path-integral kernel formula exists. Theorem 1.3 below 
provides a direct link between ensembles of non-intersecting paths and path-integral kernel 
determinant formulas, and its proof boils down to the Lindstrbm-Gessel-Viennot Lemma 
(recorded below as Lemma 2.1). By first proving the path- integral kernel determinant 
formula and then relating it to an extended kernel determinant formula, this provides 
another proof of the extended determinant point process structure for these ensembles 
(i.e., the Eynard-Mehta theorem), see Section 1.2.1. 

Let us now introduce the ensembles of non-intersecting paths. Fix T £ Z>o and N € Z>i. 
Let G = (V, E) be a finite directed acyclic planar graph with vertex set V = VoLlViU - • - UVy 
(here U represents the disjoint union of sets) and directed edge set E = Eq^\ U E\^2 U 
• • • U Ex-i^t where E n ^ n+ i only contains edges from x — > y with x € V n and y € V n +i- 
Here x — > y denotes an edge directed from x to y. 

As an example, let V n be vertices of Z 2 of the form (n, i) for n—i = mod 2 (and \i\ < M 
for some M) and let E n contain all directed edges from (x,n) to (x ± l,n + 1). Paths in 
this directed graph are trajectories of simple symmetric random walks (constrained to stay 
within distance M from the origin), cf. Figure 1. 

We define a path ir as a sequence of edges (eo = xo — > x±, e\ = x\ — > x% . . . , &r—i = 
xt-i — > xt) where Xj 6 Vj for < i < T. For such a path, let 7r(n) = x n denote the n th 
vertex in the path. 

To edges e S E we associate weights w e £ R, and to a path ir we associate a weight w(tt) 
given by the product of the weights w(e n ) along the edges e n of n. For x £ Vq and y £ Vt 
we define a transition matrix 

W(x,y) := 

ir:x—>y 

where the summation is over all paths tt from x to y. Instead of considering just a single 
path, we may consider ensembles of N non-intersecting paths from elements of Vo to 
elements of Vt (by which we mean paths which use disjoint collections of vertices). We 
define the collection of all such paths as 



Af.I.(N) := < II = {tti, . . . , 7tat} : Vi, tti goes from Vq to Vt and no two paths intersect 
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We will describe a measure on such an ensemble. This requires the introduction of two 
additional families of functions. For N fixed, consider functions ipi : Vq — > K, 1 < i < N, 
and ipj : V T -> R, 1 < j < N. Define the weight of II G M.l.(N) as 

:= det[^(Ti(0))]£-=i (n^)] det[^(^(r))]5=i, (6) 



and the partition function as 



=i 



M(n). 



ueAT.x.{N) 



If ^ and are (^-functions then the ends of the paths are fixed. 

Assuming that Z^Owe may define a measure (not necessarily positive but with total 
integral 1) on II G M.X.(N) as 

m - MS. 

When each of the three factors on the right-hand side of (6) are positive (and thus, in 
particular, v is a probability measure), one may think of v as follows. The determinants 
det[ijji(TTj(0))]fj =1 and det[(fo(7Tj(0))]^ =1 define measures on the collections of N initial 
points in Vq and N final points in Vt- The weights w{TVi) in the middle factor in (6) 
describe the transition probabilities for N independent paths (m, . . . , 7r/v) connecting these 
points. The measure v is restricted to non-intersecting paths, and the division by the 
normalizing constant Z means that v corresponds to a measure conditioned on the N 
paths not intersecting. 

Define ipf 1 : V ->■ R by <pf\x) := £ y6 y T W{x, V^jiv)- Note that this implies that W 
has a right-inverse on spanj^ ^}^ which is given by W -1 ^- = ipj. 

We will make the following biorthogonality assumption on the {ipi}fL l and {^p^}jLi- 

Yl M x )<pf\ x ) = x i=r 
xev 

Remark 1.2. Assuming Z / 0, one can show that it is always possible to perform a 
linear transformation in spanl^}^ in such a way that the biorthogonality assumption is 
satisfied and v remains unchanged. 

The final concept we introduce is that of a path-integral functional, which is any function 

/: E ^i x E^ 2 x • • • x E T -i^t ->■ R 

such that 

T-l 

/(e ,ei, . . . ,e T -i) = \\ /«(e„) 

n=0 

for functions f n : E n ^. n +i — > R. This definition extends to a directed path ir from Vb to Vt 
by setting f(n) equal to / applied to the ordered sequence of edges in n. From the function 
/ define a second set of edge weights w e := f n (e)w e were n is such that e G E n _> n+ \. With 
respect to these weights {w e }eS-E define a transition matrix 

W(x,y) = uJ(tt). 

7r:x— >y 

The following theorem is a consequence of Theorem 2.2, which is a similar result for a 
more general graph setting. The theorem shows how the path-integral kernel determinant 
naturally arises from ensembles of non-intersecting paths. A proof of the below result 
appears in Section 2. 
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Theorem 1.3. For any path-integral functional f as above 

N 

£ II /(Ti)f(n) = det(J - K + wr^jy^, 

n=(7ri,...,7r JV )eA/'.x.(jv) «=i 

where K: L 2 (Vq) — > L 2 (Vq) is given by its kernel 

N 
i=l 

As we will explain in the proof of Corollary 1.4 below, the above result can also be 
seen as a consequence of Theorem 1.1 and the known determinantal structure for ensem- 
bles of non-intersecting paths. Instead, we provide a direct and simple linear algebraic 
proof of Theorem 1.3 using only the Lindstrom-Gessel-Viennot Lemma. This provides an 
explanation of the appearance of path-integral kernel formulas. 

1.2.1. Recovering the determinantal structure. As an application of Theorems 1.1 and 1.3, 
let us see how to recover the determinantal structure of the ensemble of non-intersecting 
paths distributed according to the law v. We would like to show that for any collection of 
vertices {x\, . . . , Xk} G V, the v- measure of the set 

{n G A/".Z.(A0 : all of the visited by paths in II j 

can be written as det[K(xi, Xj)]% =1 for some fixed matrix K with rows and columns 
indexed by the set of vertices V. This property can be seen as a consequence of Corollary 
1.4 below which we show. 

Consider any collection of functions q n : V n — > R, < n < T — 1. Consider 
the space of matrices with rows and columns indexed by V, and for notational conve- 
nience denote i £ 14 as (n, x) so that matrix elements of a matrix M are written as 
M(n,x;m,y). Define a matrix Q so that Qf(n,x) = q n (x)f(n,x). For m < n and 
x £ V m , y G V n define W m;n (x, y) := J2-K-x^y w ( 7T ) (f° r m = n let this be the identity 

matrix). For x G V n define <Pj(x) := ^2 ye v T Wn,T(x,y)<Pj(y)- For x\,x 2 G V n define 
K n (xi,X2) := YliLi <Pi faltyii^) • Note that for m < n, Wm,n has a right-inverse on 
span{(^^ m ^}^ 1 which is given by W^ n ^ = <Jf^- We will write this inverse as W n , m . 
On account of this we may define the following (extended kernel) matrix 



K Qxt (m,x;n,y) 



W m>n K n (x,y) ifm>n, 
-W m ,n(I ~ K n )(x,y) ifm<n. 



Corollary 1.4. For any collection of functions q n : V n M, < n < T — 1 

N T-l 

E II II U^(n)Mn) = det(J - QK cxt ) L2(v) . (7) 

n=(7ri,...,7r A r)eA/'.X.(A r ) i=l n=0 

where we recall that ir(n) denotes the vertex in V n through which tt passes and that q(x) = 
l-q(x). 

This result is essentially a version of the Eynard-Mehta Theorem. 

Proof. We use Theorem 3.3 (the more general version of Theorem 1.1). The technical 
assumptions are immediately satisfied since we are dealing with a finite vector space. The 
right-invertibility, semigroup property and reversibility relation are all readily checked from 
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the definitions of the W mjn and K n . As a consequence of that theorem we find that we 
may rewrite the right-hand side of (7) as 

RHS(7) = det(J - K + Q Wo,iQ 1 W 1 , 2 • ■ ■ Q T - 1 W¥-i,TWr, ^o) £a(v&) , (8) 

where for x G V n , Q n f(n,x) := q n (x)f(n,x). 

Define a path-integral functional / so that for an edge e S En-+ n +i from x £ E n to 
y € E n+ i, f n (e) = q n (x). As a consequence, for any path tt from Vb to V7 1 , /(tt) = 
rin=o 9«( 7r ( n ))- This observation and Theorem 1.3 then imply that we may rewrite the 
left-hand side of (7) as 

LHS(7) = det(/ - #0 + Wo,TWT,oKo) L2{Vo y (9) 

Note that we have introduced the subscripts on the right-hand side to be consistent with 
the notation introduced before the statement of the corollary. Due to the specific type of 
path- integral functional /, it is now straight-forward to see that 

VVo,T = Q Wo,iQiWi,2 • • • Qt-iWt-i,t- 

This implies that the right-hand sides of equations (8) and (9) match and therefore com- 
pletes the proof of the corollary. □ 

1.3. Outline. In Section 2 we prove a result about ensembles of non-intersecting paths 
on directed graphs (which implies Theorem 1.3 above). In Section 3 we prove a general 
result (which implies Theorem 1.1 above) showing the equality between certain extended 
kernel and path-integral Fredholm determinants. In Section 4 we apply this equality be- 
tween Fredholm determinants to a variety of extended kernels from the literature, and in 
Appendix A we check the technical assumptions necessary in order to do this. 
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AB was partially supported by the NSF grant DMS-1056390. IC was partially supported by 
the NSF grant DMS- 1208998 as well as by the Clay Research Fellowship and by Microsoft 
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Natural Science and Engineering Research Council of Canada and by a Fields-Ontario 
Postdoctoral Fellowship, as well as by Fondecyt Grant 1120309 and Conicyt Basal-CMM 
project. DR is appreciative for MIT's hospitality during the visit in which this project was 
initiated. 

2. NON-INTERSECTING DIRECTED PATHS ON WEIGHTED GRAPHS 

2.1. A general combinatorial result. Let G = (V,E) be a finite directed acyclic planar 
graph with vertices V and edges E. Fix source vertices X = {x\, . . . , xn} and sink vertices 
Y = {yi, . . . , yjv}- Fix edge weights w e for each directed edge e G E. 

A directed path tt is a sequence of vertices connecting a vertex in X to a vertex in Y 
via directed edges in E. Denote the source vertex of ir by ir(b) and the sink vertex of tt 
by Tr(d) 1 . We say that two paths intersect if their vertex sets have non-empty intersection. 
To a directed path tt we associate a weight w{tt) which is given by the product of w e over 
edges e of tt. Define 

W(x,y) = ]T w(tt) (10) 

tt:x— >y 

where the sum is over directed paths tt from source vertex x to sink vertex y. 



We use b to denote the source, or base vertex; and d to denote the sink, or destination vertex. 
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X 7T 2 (6) 




Figure 1. An example of a graph G with source vertices X on the left- 
hand side and sink vertices y on the right-hand side. Here there are N = 3 
non- intersecting paths which are shown in grey, with starting points 7Tj(6) 
and ending points 7Tj(d) for i = 1, 2, 3. 

Define the ensemble of N directed non-intersecting paths from the elements of X to the 
elements of Y as 

J\f.l.(N;X — > Y) = {{vri, . . . ,7Tat} : Vi, 7Tj(6) £ X, 7Tj((f) E F and no two paths intersect}. 

Write n = {7ri, . . . ,7rjv} for an element of N.1.(N;X — > Y). Note that the 
non-intersection condition and the fact that X and Y have N elements ensures that 
{TT 1 (b),...,7r N (b)} = X and { 7 ri(d) ) ... J 7ri V (d)} = y. 

Lemma 2.1 ([25,29,31,44]). Fix iV > 1. For any finite directed acyclic planar graph G 
with source vertices X = {x\, . . . , xn}, sink vertices Y = {yi, . . . , yjy} and edge weights w e 
for each directed edge e G E, 

N 

det[w( Xi , yj )]^ =1 = yi ino- 

U^MX.(N;X^Y) i=l 

Consider now finite sets of source vertices X C V and sink vertices y C V (with at least 
N vertices in each of X and y). For such X and y we can likewise define the ensemble 
of iV directed non-intersecting paths from elements of X to elements of y. Denote this 
ensemble M.1.(N; X -+ y). 

Fix functions tpi : X — > M for 1 < i < N and functions ipj : y — > R for 1 < j < 
N. For II = {tti, . . . ,7T7v} with source vertices {ni(b), . . . ,7^(0)} c X and sink vertices 
{iTi(d), . . . , 7T7v(^)} C y we define the weight of II as 

Wt(U) :=det[^(&))]5-=i {j\ w ^ ^lVi(*M)]Z=v 

Define a partition function for directed non-intersecting ensembles of N paths from X to 
y with respect to weights {w e } e& E and functions {tpi}^, {(fj}^ =l as 

z = z(x,y,{w e } eeE ,{^ i }fL 1 ,{ ( p j }f =1 ):= ]T wt(u). 

IleM.X.{N;X^y) 

For 1 < j < N define ipf ] : X ->• R by 

vey 
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and further define the operator K : L 2 (X) — > L 2 {X) by its kernel (which is just an X- 
indexed matrix) 

N 

K{x 1 ,x 2 ) = Y J vt\xMi{x2)- (12) 
i=i 

Observe that W has a right-inverse on spanjt/?^}^, which we will denote by W _1 , given 
by 

In particular, since the range of K is contained in spanjt/?^}^, W~ 1 K is well defined as 
an operator mapping L 2 (X) to L 2 {y) : 

N 

w- 1 Kf = Y,(^J)m X) ^ (13) 

i=l 

where (•, •) ip.{v\ is the inner product in 1?{X ). 

We say that the biorthogonality assumption is satisfied if 



ibi, ipf ) = li-j for all 1< i, j < N. 



Theorem 2.2. Let G = (V,E) be a finite directed acyclic planar graph. Fix sets of source 
vertices X C V and sink vertices y C V. Fix edge weights w e and a second set of weights 
w e for each directed edge e G E. Fix functions tpi : X — > R for 1 < i < N and functions 
ifj : y ~R for 1 < j < N which satisfy the biorthogonality assumption with ipj defined 
via the w e weights. Write 

z = z(x,y,{ We } eeE ,{^ i }f =1 ,{i Pj }^ =1 ) and z = z(x > y,{u, e } eeE ,MiLi,{w}jLi)- 

Then 

^=det(I-K + WW- l K) L2(x) , 

where VV : L 2 {y) — > L 2 (X) is given by (10) with w replaced by w, and W^K : L 2 (X) — > 
L 2 (y) is defined in (13). 

Remark 2.3. As will be clear from the proof, the biorthogonality assumption implies that 
Z / 0. Conversely, one can show that if Z ^ then there exists a linear change of basis in 
the space spanned by {ipi}^Li ( or equally well the space spanned by {ipj}^ =l ) which leads 
back to the biorthogonality assumption being satisfied and does not change the ratio Z/Z. 

Before turning to the proof of the theorem, let us check that it implies Theorem 1.3. 

Proof of Theorem 1.3. Recall that for the path-integral functional /, we defined a set of 
weights w e = f n (e)w e where e G E n -+ n+1 . Let Z = z(v ,V T ,{'We}eeE,{' l Pi}i=vWj}f=i 

and Z = Z(Vo ) VT,{v} e } ee E,{ipi}^i,{<Pj}jLij- We claim that 



/ 



Z 



dv{Ii) f(m) = - = det(J -K + WW-^v^y 
ueAf.x.(N;V ^v T ) 1=1 
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The second equality is an immediate corollary of Theorem 2.2. To see the first equality 
above observe that 

N T-l 



z 
z 



UeAf.l.(N;V ^V T ) i=l n=0 



N T-l r N 

»r,, T . \ i=ln=0 „^„,^,{ T „ T , % i=l 



n6A/".X.(Af;yo^-V T ) 



nGA^.X.(Af;Vb^-Vr) 



as desired. 



□ 



Proof of Theorem 2.2. The proof is linear algebra. We may rewrite Z by first summing 
over the subsets of X and y which host the source and sink vertices, and then considering 
all non-intersecting paths between these sets. Thus 



Z = E det[V>j(£j)] 



N 



N 



X={xi,X2,--- ,x N }cX 
Y={y!,y2,--- ,y N }^y 



. neAT.x.(Ar ; x^y) i=i 



E det^fo)]^ (iVIdetfw^,^) 



N 

»,J=1 



det[^i(%-)]£- = i. 



The second line follows by an application of Lemma 2.1. 

We may now apply the Cauchy-Binet identity twice. The first application is with respect 
to the summation in the x's, and it yields 

iV 



v 



det 



EV'iWW^yj) 



E det[V)i(^)]-j =1 det W(xi,%) 
The second application likewise is applied to the summation in y's and yields that 



*J=1 



Z = det 



x£X 

yey 



(14) 



Observe that by the same argument we can obtain an analogous expression for Z with 
W replaced by W. However, by the definition of (p^p we find that 



Z = det 



E fatovj ( a 

.xex 



N 



«J=1 



By the biorthogonality assumption, 

,0), 



(6) 



>-l=J, 



J /L 2 (A-) 

and hence Z = 1. Now observe that we can rewrite the kernel in appearing in (14) as 



x£X 

yey 



by using the definition of the operator W and the inner product on L 2 (X). 
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det(/ -K + WW~ l K) L 2 



(X) 



det 



N 



(15) 



To prove the above statement we will write down the matrix for the operator I — K + 
WW~ 1 K in the basis (tpi , • • • , tp^ , (spanl^}^!) -1 ), where (spanltpi}^) 1 - represents any 
basis of the orthogonal complement of spanj^j}^ in L 2 (X). Let us consider the action 

of the operator K on the basis elements. On tpj one sees that K acts as the identity 
operator: 

x 2 £X \i=l / i=l \x 2 eX 



i=l 



(6) 

J / L 2 (X) 



N 



(6), 



i=l 



It is likewise clear that K acts on the basis elements of (spanl^}^) -1 - by taking them all 
to zero. Thus we may write K as the matrix 



K 















where the two blocks correspond to the basis elements {ip^ }jLi an d (span{^j}^ 1 )- L . This 
shows that 

I — K 















The remaining operator to study is WW 1 K. Writing the corresponding matrix in 
blocks as above we get from (13) that 



A 





* 






where the N x N matrix A is yet to be determined. The value of the star is not important. 
To see this, write 



i- k + ww -1 ir 



A 





* 





and observe then that 



det (J -K + WW~ l K) L 2 [x) = det^] £ =1 



(16) 



The value of can be found by using the inner product, 



Recalling that W 1 Kip^ = <pt we deduce that 



.4 



i.j 



L\X) 



Combining this with (16) proves (15) and hence completes the proof of the theorem. □ 
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3. Equivalence of extended kernel and path-integral kernel Fredholm 

determinants 

There are various types of limits one can take of graph-based non-intersecting line en- 
sembles. In this section we will show that formulas of the type given in the previous section 
survive these limits. We do not prove this directly via a limit transition, but rather show 
how such formulas arise via manipulations of the extended kernel Fredholm determinants 
which describe these limiting systems. The main result of this section is, therefore, the 
equality of two types of Fredholm determinants. 

This equality will be stated in an abstract setting in this section, and later applied to 
the examples we are interested in in Section 4. A concrete example to keep in mind is the 
Airy2 process, which we will use throughout this section to illustrate the objects we will 
introduce and the assumptions we will make on them. The Airy2 process was introduced 
in the Introduction and is discussed in further detail in Section 4.2, we refer the reader 
there for details and just recall the definitions of the Airy Hamiltonian H = — A + x and 
the Airy kernel K-zix, y) = J Q dX Ai(x + A) Ai(y + A). 

3.1. Fredholm determinants. Let us briefly introduce some of the basic notions related 
to Fredholm determinants (we refer the reader to [43] for more details). Consider a separa- 
ble Hilbert space H and let A be a bounded linear operator acting on % (H = L 2 (R) in the 
Airy 2 case). Let \A\ = V A* A be the unique positive square root of the operator A* A. The 
trace norm of A is defined as ||A||i = Y^=i( e m \A\ e n), where {e n } n >i is any orthonormal 
basis of %. We say that A G B\(J-L), the family of trace class operators, if ||^4||i < oo. 
For A G £>i(%), one can define the trace tv(A) = ^2^ = i{e n , Ae n ) . For later use we also 
define the Hilbert- Schmidt norm \\A\\2 = \/tic(\A\ 2 ) and say that A G B^^H), the family 
Hilbert- Schmidt operators, if ||^4||2 < oo. Given A G Bi(H) one can define a generalization 
of the finite-dimensional determinant, the Fredholm determinant det(/ + ^4)%. We refer the 
reader to [43] for the details of the definition in this level of generality and just point out 
that, as expected, det(I + A)% = \\ n {l + A n ), where A n are the eigenvalues of A (counted 
with algebraic multiplicity). 

The result presented in this section (Theorem 3.3) can be stated, under some conditions, 
for operators acting on a general separable Hilbert space. Nevertheless, in order to keep 
the presentation as simple as possible, and since it is the setting we need for the examples 
in Section 4, we will restrict ourselves to the case of integral operators on an L 2 space. 

More precisely, we assume we are given a measure space (X, X,/z) and consider the 
Hilbert space L 2 (X,fi). For brevity we will drop // from the notation. We will also denote 
by Ai(X) the space of real- valued measurable functions on X. By an integral operator 
we mean an operator A: D C M.(X) — > M(X) acting as Af(x) = f x n{dy) A(x,y)f(y), 
where A: X x X — > R is the integral kernel of A. We will often speak interchangeably 
of an integral operator and its kernel. In particular we have abused notation by using the 
same letter to denote an integral operator and its kernel. We recall that the product of 
two integral operators is defined by AB(x, y) = j x fi(dz) A(x, z)B{z, y). 

Though we will not appeal to this, we note that the Fredholm determinant det(/ — 
K) L 2(x) °f a trace class operator K : L 2 {X) — > L 2 (X) with continuous (in both x and y) 
integral kernel K{x,y) has the following (absolutely convergent) series expansion 

det(J - K) V 2 (X) = l + J2 I d ^ • • ' / Mxk) det[K( Xi , Xj)$ j=1 . (17) 

fc>i Jx Jx 

3.2. Assumptions for the theorem. In order to state the main theorem of this section in 
a fairly broad context, we must introduce a few operators and impose certain assumptions 
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upon them. Most of the assumptions are technical and intended to ensure well-definedness 
or finiteness of the various quantities involved in the statement of the theorem. The main 
(not just technical) assumption is given in Assumption 2. 

Fix t\ < • • • < t n for the duration of this section. We will be interested in comparing 
the Fredholm determinant of certain integral operators acting on the Hilbert spaces L 2 (X) 
and L 2 ({ti,...,t n }xX) (the measure we use in the second space is the product of the 
counting measure on {t\, . . . ,t n } and //). The operators we consider will be constructed 
from the following four families of operators: 

• For each 1 < % < j < n, an integral kernel Wt^tj (for convenience we also introduce 
the notation = /); 

• For each 1 < i < n, an integral kernel K^; 

• For each 1 < i < j < n, an integral kernel W^^K^ (for convenience we also introduce 
the notation W tai K ti = K u ); 

• For each 1 < i < n, a multiplication operator Q% i acting on Ai{X) as Qtif{x) = 
q u (x)f(x) for some q u G M(X). 

The reason for the choice of notation Wt in the third family of operators is that we 

will assume below that W^^W^.^iftj = for i < j (so that even though it is not defined 
above as its own operator, VVV t- can be thought of as a right inverse of WV f ° n the range 
oiK ti ). 

We make the following (technical) assumption. 
Assumption 1. 

(i) The integral operators Q ti W ti , % , Q ti K u , (^.VV/ ./ A'/ ; and Q tj W tj ,uK u for 1 < i < 
j < n are all bounded operators mapping L 2 (X) to itself. 

(ii) The operator 

K tl - Q tl W tl ,t 2 Q t2 ■ ■ ■ Wt„_ 1)tn g tn W tn)tl if tl , 
where Q t . = I — Q^, is a bounded operator mapping L 2 (X) to itself. 

The last operator in the assumption will appear in the formula provided in Theorem 
3.3. An alternative expression for this operator, which is in some cases more convenient 
for checking the assumption, is given in Lemma 3.1. 

In the case of the Airy2 process we take X = M, choose /i to be the Lebesgue measure 
and set W ti , tj = e^~^ H , K u = K 2 , and W tj , h K u = e^~^ H K 2 for 1 < % < j < n. 
One can take for example the operators to be projections on intervals [a>i, oo), that is, 
Quf( x ) = ^-x>aif(x), which corresponds to studying the finite dimensional distributions 
of the Airy2 process (we will make a more general choice in Section 4.2). 

Going back to the general setting, we will make a certain algebraic assumption on the 
operators V\? ti ,t v K u and W tht .K ti . 

Assumption 2. For each i < j < k the following hold: 

(i) Right-invertibility: yV t .j.yV t . tt .K ti = K ti ; 

(ii) Semigroup property: W ti , 4j W 4j ,t fe = W ti ,t fc ; 

(iii) Reversibility relation: j . A'/ = K u W tl ,t r 

The second property is clear in the Airy2 case, while (i) and (iii) follow from the fact 
that K 2 is the projection operator into the negative (generalized) eigenspace of the Airy 
Hamiltonian H (see Section 4.2). 
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Let us now explain how these operators will be used. Using the kernels introduced above 
we define an extended kernel K cxt as follows: for 1 < i,j < n and x,y E X, 

(W t t Kt (x,y) if z > ?, 

K^(t i ,x;t j ,y) = \ ^ ); ~ J (18) 

[ -Wu,t 3 K 1 ~ K t 3 ){x,y) if * < j . 

This definition coincides with the usual notion of extended correlation kernels of determi- 
nantal point processes, cf. [7,17,23,26,33,45]. In the case of the Airy 2 process, it coincides 
with the definition given in the Introduction and in (40). As an operator, K cxt acts on 
/ e Ll c ({h,...,t n }xX) as 

n „ 

3=1 Jx 

See Section 4 for concrete examples. 

We also need to make the following (technical) analytical assumption. 

Assumption 3. One can choose multiplication operators V^, Vp, Ut t and Up acting on 
M(X), for 1 < i < n, in such a way that: 

(i) VlVuQu = Qu and K ti U' t JJ ti = Kp, for all 1 < i < n. 

(ii) Theoperators V ti Q ti K ti V{ v V u Q u W u ,t 3 V( p V u Q u W u , tj K t] V^ and^.W^i^ 
preserve L 2 (X) and are trace class in L 2 (X), for all 1 < i < j < n. 

(iii) The operator U ti [w^K^ - Q ti Wt it t i+1 ■ ■ ■ Qt^t^nQtW^Kt,]^ pre- 
serves L 2 (X) and is trace class in L 2 (X), for all 1 < i < n, where Q ti = I — Qp. 

The primes in U[. and V(. mean that these are almost (left) inverses of the operators Up 
and Vp , and hence the multiplication by these operators in (ii) and (iii) should be thought 
of as a conjugation. The distinction is because in many cases it will be necessary to let Vp 
be multiplication by a function which is where qp is 0, in which case Vp is not invertible, 
with an analogous situation for Up and Up. 

Before stating the main result of this section, Theorem 3.3, let us state a formula which 
reexpresses the operator appearing in Assumption 3 (iii). Besides being used in the proof 
of the below theorem, this formula is often useful in checking the assumption (for example, 
as in Remark 3.2). 

Lemma 3.1. Writing Q t = I — Qt, we have, for any 1 < i < n, 

Wp M K tl - QpWp, t2 Q t2 ■ ■ ■ Wtn-utnGtnWtnfrKt! 
n n—j 

j=i k=0 j=ao<ai<---<a k <n 

We postpone the proof of this lemma until the end of this section. 

Remark 3.2. Suppose that there exist multiplication operators Vp and Vp acting on 
M(X), for 1 < i < n, in such a way that: 

(i) KVpQp = Qp and KpVpV^ = K t ., for all 1 < i < n; 

(ii) Theoperators V ti Q ti K u Vl v VpQpW u , t] V{ p V u Q u W u , tj K tj V{. and % Q tj W tj , k K u % 
preserve L 2 (X) and are trace class in L 2 (X), for all 1 < i < j < n. 

Then it is not hard to check, using the formula given in Lemma 3.1, that Assumption 3 
holds, taking Up = Vp = V{. and Up = Vp = Vp (see the end of the proof of Corollary 4.6 
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in Appendix A for more details). In the case of the Airy2 process, when the operators Qt t 
are of the form Q ti f(x) = l x > a .f(x) as discussed above, both V ti and V{. can be taken to 
be the identity. If, on the other hand, one assumes ^(x) to be for x < but to grow 
at a certain rate for x > ai, as we will in Section 4.2, then it is necessary to choose these 
operators more carefully (see the proof of Corollary 4.6). 

3.3. Identity between extended and path-integral kernel Fredholm determi- 
nants. Define a diagonal operator Q acting on / G A^({ii, . . . ,t n }xX) as 

Qf(U,-) = Quf(ti,-)- (19) 

Note that, by Assumption 1, QK ext preserves L 2 ({ti, . . . , t n } x X). The following result 
expresses the Fredholm determinant of / — QK cxt on L 2 ({ti, . . . ,t n } xX) as a Fredholm 
determinant on L 2 (X). The first example of such a formula was provided by [36] for the 
case of the Airy2 process (see also [37]). This was later extended to the Airy! process in 
[40]. This type of formulas have recently been found to be very useful in the study of these 
processes, see for example [22,32,38-40]. 

Theorem 3.3. With the above notation, and under Assumptions 1, 2 and 3, we have 
det{l-QK^) L2[{ti M ^ 

(20) 

where Q t = I — Q t . 

Remark 3.4. The operators appearing in both Fredholm determinants preserve L 2 (X) by 
Assumption 1. Moreover, the Fredholm determinants are well-defined thanks to Assump- 
tion 3, even though the operators appearing there are not necessarily trace class. In fact, 
if we define the diagonal operator V acting on u G L 2 {{t\, . . . , t n } xX) as (Vu)^ = V^u^, 
and similarly define V , then VQK ext V' is trace class by Assumption 3(h) and by the cyclic 
property of the determinant and the fact that V'VQ = Q it leads to the same Fredholm 
expansion for det(J- VQK eid V') L 2 i{tli ... >tn}xX) and for det(/-QK cxt ) L 2 ({tlv .. jin}xX) . The 
same argument applies to the Fredholm determinant on the right-hand side of (20) by 
Assumption 3(iii), if we multiply it on the left by Ut l and on the right by U[ . Hence both 
sides of (20) are well-defined and one should really read the equality as 

det(l- VQK ext V') r2fU . , v , 

V / L 2 ({ti,...,t n }xX) 

= det(J - U tl (K tl - Q tl W tl , t2 ■ ■ ■ Wt n - 1 ,tnQ tn Wtn,t 1 KM 1 ) I ? {x) . 

Proof of Theorem 3. 3. The proof of this result is a generalization of the proof of Theorem 
1 of [40] (see also the Appendix of [37]). We will retain most of the notation of [37,40], and 
as in those papers we use sans-serif fonts (e.g. W) for operators on L 2 ({t\, . . . ,t n } xX). 
This space can be identified with the space © tg { tl tn yL 2 (X), and hence we may (and 
will) think of an operator W on L 2 ({t\ , . . . , t n }xX) as an operator-valued nxn matrix. We 
will use serif fonts for the matrix entries (e.g. Wjj = W for some W acting on L 2 (X)). All 
determinants throughout this proof are computed on L 2 ({ti, . . . ,t n }xX) unless otherwise 
indicated. 

We will use repeatedly the following facts about trace class operators and Fredholm 
determinants on a separable Hilbert space %: 

(i) If A, B G Bi(H) then AB € Bi(H) and 

det((J + A)(I + B)) n = det(J + A) H det(J + B) n . 

Moreover, if A and B are bounded linear operators on H and both AB, BA G Bi(H) 
then 

det(I + AB) H =det(I + BA) H . (21) 
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(ii) An operator acting on ® t6 j 4l t \ H is trace class if and only if all of its matrix 
entries are trace class. 

To simplify notation throughout the proof we will replace subscripts of the form ij by i, so 
for example Wjj = W tiitj . 

Recall that we are assuming t\ <t2 < ■ ■ ■ <t n . Let K = QK cxt . Then K can be written 

as 

K = Q(V\TK d + W+(K d -l)), (22) 

where 

and W _ , W + are lower triangular, respectively strictly upper triangular, and defined by 

w v n-,,i, v , \N± = w id i i<r 

Here we are slightly abusing notation, because Wjj is not defined for i > j. However, since 
W~ appears applied after K d , the formula makes sense, with [W-K d ] M - = WijKj for i > j. 
We also define the diagonal operators V, V, U and U' by 

Vjj = Vil i=j , V-j = V[l i=j , Uij = Uil i=j and ll'^ = U-l i=j . 

In order to manipulate the Fredholm determinant of I — K we will need to make sure at 
each step that the appropriate operators preserve L 2 {X) and are trace class in L 2 {X) as 
needed. As a consequence, the proof is slightly cumbersome, so we will first briefly explain 
the main idea, ignoring some details and all analytical issues. 

Our goal is to manipulate the determinant of I — K in such a way that we end up with 
the determinant of an operator-valued matrix I — K where only the first column of K is 
non-zero. If we achieve this, then we will have det(l — K) = det(l — K) = det(l — ^x,i)l' 2 (X)i 
and all we will need to do is compute Kn. The key to obtain such an identity is the 
following observation. Using the semigroup property in Assumption 2(h) one can check 
directly that 

[(I + W+)" 1 )]. . = Il i= i - W i)i+ il,- =i +i. (23) 

This identity is meant in the sense of products of integral kernels, where the product of 
the identity operator with an integral kernel is defined in the obvious way. Now using the 
identity Wij-iKj^iWj-ij = WijKj from Assumptions 2(h) and 2(iii) we get that 

[( W - + W+)K d (l + W+)" 1 ] . . = WijKj - Wij-iKj-iWj-ijlfri = Wi,iKxl j=1 . (24) 

Note that only the first column of this matrix has non-zero entries. To take advantage of 
this fact we rewrite K as 

K = Q(W~ + W+)K d (l + W+)- 1 (l + W+) - QW+, (25) 

so that 

I - K = (I + QW+) [I - (I + QW+)^ 1 Q(W^ + W+)K d (l + \N + y\\ + W+)] . 

The invertibility of I + QW + follows from the fact that QW + is strictly upper triangular. 
This fact also implies that det(l + QW + ) = 1, and hence 

det(l - K) = det(l - (I + QW+)- 1 Q(W- + W+)K d (l + W+)- 1 (l + W+)) 

= det(l - (I + W+)(l + QW + )- 1 Q(W~ + W+)K d (l + W+)" 1 ), 

where we have used the cyclic property of the determinant. Recalling that only the first 
column of (W~ + W + )K d (l + W + ) _1 is non-zero we deduce that 

K = (I + W+)(l + QW+)^ 1 Q(W" + W+)K d (l + W+)" 1 

has the same property and hence det(l — K) = det(l — K) = det(l — Ki^)^^) as desired. 
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The rest of the proof will consist in making the above argument rigorous and precise 
and then computing the resulting Ki i. Recall that, by Assumption 3(h), each entry in the 
operator- valued matrix VKV is trace class in L 2 (X). Let 

Wi = VQW+V and W 2 = VQ(V\T + W+)K d V. (26) 

Since VQW+V" is strictly upper triangular, we have (VQW+V) n+1 = 0, so I + Wi is 
invertible: 

n 

(I + Wi)" 1 = ^(-l) fc (VQW + V / ) fc . (27) 

k=0 

Therefore we can write 

det(l — VKV") = det((l + Wi)(l - (I + Wi)~ 1 W 2 )). 

We remark that Wi, W 2 and (I + Wi) -1 are trace class in L 2 (X) by Assumption 3(h) and 
(27), and thus from the last identity we deduce that 

det(l-VKV') = det(l +Wi) det(l - (I + Wi)- 1 W 2 ) = det(l - (I + Wi) _1 W 2 ) , (28) 

where the second equality follows from the fact that, since Wi is strictly upper triangular, 
its only eigenvalue is 0, so det(l + Wi) = 1. 

Write 

(l+Wi)- 1 W2 = W 3 W 4 (29) 

with 

W 3 = (I +Wi)- 1 VQ(W- + W+)K d (l +W+)- 1 U / and W 4 = U(l + W+)V. 

Here we are using (24) and the identity KUU' = K. We have already checked that W3W4 is 
trace class in L 2 (X). Thus if we prove that W4W3 is also trace class we can deduce from 
(21), (28) and (29) that 

det(l-VKV') = det(l — W4W3). (30) 

We want to obtain an explicit expression for the kernel W4W3. Note that, in view of 
(27) and the fact that V'VQ = Q, V'(l + Wi)^VQ = (I + QW+)^Q, so all the factors V 
and V cancel in \N4\N3: 

W4W3 = U(l + W+)(l + QW+)~ 1 Q(W~ +W+)K d (l + W+^U". 

Prom (24) and the semigroup property we deduce that, for < k < n — i, 

(QW+) fc Q(V\T + W+)K d (l + W+)" 1 

J i,l 



i<oi<' ■ ■ <at<n 



while for k > n — i the left-hand side above equals (the case k = is interpreted as 
QiWi^iKi). Summing the above times (— l) k from k = to k = n — i we get directly from 
the last formula and (27) that 



(I + QW + )" 1 Q(W" + W+)K d (l + W+) 



1 j 



i;=i 



QiW^K! + ^("l) fc QiW^lQa^a^ ■ ■ ■ Qa k _ L Wa k ^a k Qa k Wa k ,lKl 

k=l i<ai<---<a k <n 



(31) 
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Note that only the first column of the above matrix contains non-zero entries. Since 
U(l+W + ) is upper triangular and U' is diagonal, the same is true for W4W3. Pre- multiplying 
(31) by U(l +W+) we get 

n n-j 

(W 4 W 3 )- 1 = ^^(-l) fc Yl UiWijQjWj^Wa^Qa^Wa^Q^W^KtUi. 
j=i k=0 j—ao<ai<--<a^<n 

(32) 

By Lemma 3.1 we deduce that 

(W4W3) . 1 = Ui [Wi.iifi - QiW i>2 Q 2 ■ ■ ■ W n -i, n Q n W n ,iXi] U[. (33) 

By Assumption 3(iii) this operator is trace class, which provides the needed justification 
for writing (30), and then since only the first column of W4W3 is non-zero we deduce that 



det(l - VKV') = det(J - (W 4 W 3 ) 1 ly 

Setting i = 1 in (33) yields the result. □ 

In order to finish the proof of Theorem 3.3 it remains to prove Lemma 3.1. 

Proof of Lemma 3. 1 . We start with the right-hand side of the identity. Replace each Q, L 
by I — Qi except for the first one to get 

n n—j k f _ ■ _ \ 

EE("!) fc E VmJ ( " ir £ W i>6o Q 6o W 6o>6l ■■■Q bm _ 1 W^^Q^W^K, 

j=i k=0 m=0 ^ ' j=b Q <bi<---<b m <n 

where, as in the above proof, we have written i instead of ti in the subscripts. Interchanging 
the order of summation leads to ^" =i S^„(-l)' !+m ( , '^~ m ) (*)j where (*) repre- 

sents the last sum above, and is independent of k. Noting that Y^k=m { n ~k-~rn l )(~^) h+m = 
l m=n _j, the above expression can be rewritten as 

n 

E E Wi,boQb m AQ bl ■ ■ ■ m^K-^K-^K-^Kx 

j=i j=b <b 1 <---<b n - j <n 
n 

= E w ^ - QiWjj+iQi+i ■ ■ ■ m^QbWn,!^ 

j=i 

n 

= E [ W ij+lQj+lWj+ld+2Qj+2 ■ ■ ■ Wbn-ubnGbn'Wn.lKl 

j=i 

- WijQjWjj+iQj+i • • • w 6n _ li6n Q 6n w re ,i#i 

= W itl K! - QiWi^Qi • • • W n _i, n Q n W n ,iKi, 
where the last equality follows by telescoping. □ 



4. A FEW EXAMPLES 



We will now show how to apply Theorem 3.3 to a few examples of Fredholm determinants 
which arise in describing objects of interest in random matrix theory, growth processes, 
particle systems, tilings and representation theory. Our examples include extended deter- 
minantal point processes such as the stationary (GUE) Dyson Brownian motion, the Airy2 
process, and the Pearcey process; all of which are limits of ensembles of non-intersecting 
directed paths on weights graphs. We also include an extended determinant point process 
given by Markov processes on partitions related to the z-measures; this ensemble is not a 
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limit of a graph-based ensemble of non-intersecting directed paths. We also show how the 
identity applies to signed extended determinantal point processes such as the Airyi and 
Airy2_i>i processes. 

The proofs of the results in this section are postponed to the appendix. 



4.1. Stationary (GUE) Dyson Brownian motion. Consider the eigenvalues of an 
NxN Hermitian matrix with each (algebraically independent) entry diffusing according to 
a stationary Ornstein-Uhlenbeck process (real valued on the diagonal and complex valued 
off the diagonal). The eigenvalues of this process are real valued and themselves form a 
Markov process, called the stationary Dyson Brownian motion. Its stationary marginal 
distribution is the NxN Gaussian Unitary Ensemble (GUE) eigenvalue distribution. 
We will consider this Dyson Brownian motion process in stationarity and write the i th 
largest eigenvalue at time t as Xi^(i;t). The collection of eigenvalues at time t is written 
Ajv( - ;i) = (Ajv(1; t), . . . , Xn(N; t)) and the curve traced out by the i th eigenvalue over 
time is written as Atv(z; •). Then the graphs of Ajv(1; •))•••! ^n(N; •) form an ensemble of 
non-intersecting curves (see, for example, Section 4.3.1 of [4]). This ensemble of curves 
is indexed by time t and curve label i and hence can be thought of as a random variable 
taking values in the space of continuous curves from {1, . . . , N} x M to M. We will write E 
as the expectation operator for this random variable. 



Definition 4.1. For times ti < t2 < 
qti(x) = 1 — qti(x)- For a curve g : E — > I 
One likewise defines the functional q(g) 



■ • < t n consider functions : 
define the functional q by q{g) 

-Y\U<nMu))- 



n 



R and let 

-.iquigiU)). 



The stationary (GUE) Dyson Brownian motion is an extended determinantal point pro- 
cess. In particular this means that for any functions q^ (as above), 



E 



N 



det(l - Q#GUE,jv)L2({t 1 ,...,t„}x] 



(34) 



as long as both sides are well-defined, where Q is defined as in (19) and -Kque n 1S ^ ne 
extended Hermite kernel (see e.g. [45]): 



K^ N (s,x;t,y) = { 



N-l 



k=0 



- eHs t] Mx)^k(y) if a <t. 

k k=N 

Here <fk(x) = e~ x2 / 2 pk(x) and is the fc-th normalized Hermite polynomial (so that 

ll^fclb = !)• 

Writing 

D = -|(A- x 2 + 1), 
the harmonic oscillator functions (fk satisfy Dpk = kipk- Then the Hermite kernel 

N-l 

Kgve,n{x,v) = K GVEjN (0,x;0,y) = ^ Pk(x)(p k (y) 

k=0 
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acts as the projection operator onto span{c/?o, • • • , yjv-i}- In the notation of Theorem 3.3 
we are taking X = R, [h the Lebesgue measure, and for 1 < i < j < n 



W u , tj (x,y) = e-^^ D (x,y) = ]T e-^-^ k y k {x)tp k {y), K u = K GVE 



N • 



k=0 



N-l 



W tj , ti K ti (x,y) = e&-^ D K GVE , N (x,y) = E e {t ^ u)k ^ k {x)^ k {y). 



k=0 



Applying Theorem 3.3 we conclude: 



Corollary 4.2. Fix t\ < ■ ■ ■ < t n and write r = min^i^.^-i — ti\. For each 1 < 
i < n choose a function £ L 1 1 oc (M) satisfying sup xeR e~ Kx2 \qt i (x)\ < oo for some k € 
(0,^tanh(r/v^)). Th 



E 



en 



u 



det(7 - Kgue.tv + Q tl e {tl - t2)D Q t2 ■ ■ ■ Q t y tn ~ tl>u K GVE 



{t„-ti)D, 



N, 



Note that we have removed the bars over q and Q by replacing g^'s by (1 — qtj's. 



4.1.1. Continuum statistics. We may now take a continuous time limit of the above formula 
(in the style of [22]). Consider a function /i:lxl->[0, oo] and I < r. Define an operator 
r^ - acting on L 2 (R) as follows: r f{-) = u(r, •), where u(r, •) is the solution at time r of 



d t u 



-Du — hu 



(35) 



with initial data u(£,x) = f(x). By the Feynman-Kac formula we may also express the 
action of this operator in terms of a path-integral through a potential h as 



(36) 



with 



where the expectation is over a (standard) Brownian motion &(•) started at time 
b(£) = x and run until time r. 

Let t\ = £, t n = r and the U be spaced equally in between with step size 5 = (r—£)/(n—l). 
Then letting qti(x) = 1 — 5h(tt,x) and taking n — > oo the above formula yields: 

Proposition 4.3. For any interval [£, r] and continuous bounded function /i:lxl-> 
[0,oo] 



E 



N 

Jexp 



h(t,\ N (j;t))dt 



det{l-K GVEtN + T^ r e (t --^ D K GVEtN ) L2 (37) 



Remark 4.4. The condition on h is not optimal, but it makes the arguments simpler. 
A different class of functions h for which the result holds is the following. Fix a function 
g G 77 1 ([£, r]) and set h(t, x) = for x < g(t) and infinity otherwise. Then the left-hand side 

of (37) becomes P D^Lii-V/vO;*) < 9{t)Vt £ {^ r W an d the right-hand side makes perfect 
sense as well, with T^ r now being the solution operator of a certain boundary operator 
involving g. This case corresponds to calculating the probability that on the entire interval 
[£,r], the top curve of the Dyson Brownian motion remains below the function g(t). This 
is the same type of result shown in [22] for the Airy2 process, and the proof for this case 
can be easily adapted from the arguments in that paper. 
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4.1.2. Rescaled process. Now introduce the rescaled process 

X N (i;t) = V2N 1/6 (\ N (i;N~ 1/3 t)-V2N). 

Changing variables x i— > ~^^TJE X + V / 2iV, y i— > ^^ 1/6 y + \/2iV in the kernel accordingly, 
we immediately obtain: 



Corollary 4.5. For any t\ < ■ ■ • < t n and functions : 
same conditions as in Corollary ^.2, we have 



1 <i <n, satisfying the 



E 



N 



det(/ - i^GUE,iv + Q tl e {tl ~ t2)HN Q t2 ■ ■ ■ Qt n ^ tn - tl)HN K GUE>N ) 



(38) 



where the kernel of Kgve,n is given by 

KGVE,N(x,y) 
and the operator 



K G ve,n ( * r . + V^V, 7=~~r, 



.1- 



The above rescaling corresponds to focusing in on the top curves of the Dyson Brownian 
motion. In the limit N goes to infinity, Kque N converges to the Airy2 kernel K2 and Hjy 
converges to the Airy Hamiltonian H (defined in the Introduction and below in Section 
4.2). So in the limit as iV goes to infinity we recover the formula for the Airy2 process 
as expected. The operator in the Fredholm determinant in the right-hand side of (38) 
converges in trace class to the corresponding one with K2 and H, which means that all 
of the left-hand side probabilities have limits. This can certainly be proved under some 
additional (though not optimal) assumptions on the as in Corollary 4.6, but we choose 
to treat the Airy2 process independently. 



4.2. The Airy2 line ensemble. The multi-layer Airy2 process [26, 36] is the limit of 
the stationary (GUE) Dyson Brownian motion under the scaling of Section 4.1.2. In 

particular for f 6 1 consider the point process corresponding to 

As iV goes to infinity, this point process converges in the vague topology to a limiting point 
process with an infinite number of simple points which we write as {Airy 2 («; t) : i € Z>i} 
(labeled so that Airy 2 (i;t) > Airy 2 (j';t) for i < j). This convergence can be strengthened 
so that for any fixed set t\ < ti < ■ ■ ■ < t n , the n-tuple of A-point processes has a 
limit {Airy 2 (i;t) : i E Z>i,i £ {t±, . . . ,t n }}. This limiting collection of point processes is 
consistent and can be completed to a point process valued stochastic process indexed by 
t G R. This process is called the multi-layer Airy2 process. As it is the limit of a stationary 
(in t) process, it is also stationary. 

There exists a continuous version of this process [21] so that Airy 2 can be thought of as a 
random variable taking values in the space of Z>i indexed, continuous and non- intersecting 
curves from R to R. The convention is that Airy 2 (l; •) represents the top curve (i.e., the 
limit of A/v(l; •)). The continuous version of the multi-layer Airy2 process is called the 
Airy 2 line ensemble. 

Since the Dyson Brownian motion was an extended determinantal point process (34), so 
too is the multi-layer Airy2 process. Analogous to (34), and with the functional q given in 
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Definition 4.1 and operator Q given in (19), 



E 



JJg(Airy 2 (i;-)) 



det(/-QK 2 cxt ) 



L^({t 1 ,...,t n }xl 



(39) 



where AT 2 xt is the extended Airy 2 kernel 



K^(s,x;t,y) = { 



/ d\ e~ A ( s ~*) Ai(x + A) Ai(y + A) if s > t, 
Jo 

r° 

- dX e ~ A ( s ~*) Ai(x + A) Ai(y + A) if s < t, 

J — oo 



(40) 



and Ai(-) is the Airy function. In order for the above expectation to make sense, one has 
to impose conditions on the functions q^, such as in Corollary 4.6. 

To put this example in the setting of Theorem 3.3 we take X = R, n the Lebesgue 
measure, and consider the Airy Hamiltonian defined as 

H = -A + x. 

H has the shifted Airy functions Ai^(a;) = Ai{x — X) as its generalized eigenfunctions: 
HA\\(x) = AAi^(x). Define the Airy 2 kernel K 2 as the projection of H onto its negative 
generalized eigenspace: 



K 2 {x,y) 



poo 

/ dX Ai(x + A) Ai(y + A). 
Jo 



Then it is not hard to check that, in the notation of Theorem 3.3, (40) corresponds to 
taking, for 1 < i < j < n, 

/oo 
dXe x ^-^M(x + X)M(y + X), K u = K 2 , 
-00 

pOO 

W t] , u K u (x,y) = e^-^ H K 2 (x,y) = / dX e ~ A ^ _ti ) Ai(x + A) Ai(y + A). 



Note that Wti,tj is only well-defined on the range of K 2 . Applying Theorem 3.3 allows to 
conclude: 

Corollary 4.6. Fix t\ < ■ ■ ■ < t n and let r = minj = i ) ... )n _i{|tj+i — t{\}. Choose functions 
Qu ^ -^ioc(^)' ^ — * — n > suc h ^at sup a;>0 e~ rx \qt i (x)\ < 00 for some < r < r and 
sup x<0 f(x)\l — qt t (x)\ < 00 for some function <p(x) such that J° dx e -2 ^™ - * 1 ^ tp(x)~ 2 < 
00, 1 < % < n. Then 



E 



JJg(Airy 2 (j;-)) 



det(/ - K 2 + Q tl e^-^ H Q t2 ■ ■ ■ Qt n e^~ tl)H K 2 ) 



L 2 ( 



This formula is also the limit of the right-hand side of (38) as N goes to infinity. 

Since the Airy line ensemble is a continuous version of the multi-layer Airy 2 process, 
we may take a continuum limit of the above formula, in the same manner as done in 
Section 4.1.1. The PDE which T^ r is solving is now dtu = —Hu — hu (corresponding to 
replacing D by H in (35)) and the result is that for any interval [£, r] and suitable function 
h : R x R — )■ [0, 00] (for example h can be taken to be bounded, continuous, and such that 
h(t, x) = for any t £ [£, r] and x < M for some M € R), 



E 



h(t,Ahy 2 (j;t))dt 



det(l-K 2 + F f l r e^ H K 2 ) 



(41) 
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We will omit the proof of this statement, which can be adapted from the proofs of Proposi- 
tion 4.3 and Corollary 4.6 together with the proof of Proposition 3.2 in [22]. Taking h(t, x) 
to be for x < g(t) and infinity otherwise we recover Theorem 2 of [22]. 



4.3. The Pearcey process. There are many other multi-layer processes which arise as 
scaling limits of non-intersecting ensembles of Brownian motions (or similar diffusions) for 
which we can apply Theorem 3.3 (see for instance Airy-like processes [2,3,6,16]; bulk limits 
such as the Sine process [45], Pearcey process [46] or Tacnode process [8]; hard edge limits 
like the Bessel process [30]). 

To illustrate this point we will show how a Fredholm determinant involving the Pearcey 
kernel can be rewritten via Theorem 3.3. 

Let us briefly and informally recall one way the Pearcey process arises as a scaling limit 
of Brownian bridges. Consider 2N Brownian bridges on the time interval [— N, N] such 
that all 2N of them start at height and N of them end at height b and the other iV end at 
height —b. Condition these Brownian bridges not to intersect (as can be done by spacing 
their starting and ending points by e and letting e go to zero). When b = the limit shape 
of the ensemble of conditioned Brownian bridges has a limit shape which is elliptical (and 
the ensemble is sometimes called a watermelon) and the fluctuations around the top of this 
limit shape are described (in the limit as N goes to infinity) by the Airy2 line ensemble 
minus a parabolic shift. 

When the endpoints parameter b = cN, the limit shape has a cusp at some time t = c'N, 
where d G (—1,1) is a function of c. For t± < ti < ■■■ < t n , the iV-tuple of point 
processes formed by the heights (properly centered and normalized by iV 1 / 4 near the height 
of the cusp) of the Brownian bridges at times c'N + tjA^ 1 / 2 , 1 < i < n, converges in the 
vague topology as N goes to infinity to a limit which is called the Pearcey process, V, see 
[5,19,20,34,46]. It is a point process valued stochastic process indexed by t S R. At each 
time t the point process can be indexed by Z as {V(j; t) : j £ Z}. 

Analogously to (34), and with the functional q given in Definition 4.1 and operator Q 
given in (19), 



E 



II ^(i; 



det(l-QK&) (42) 



j=-oo 

where K^f c is the extended Pearcey kernel 

+ Wfj c dU j_,J V e-^/^,2- X u v _ u ' ( 43 ) 

and where C is the contour consisting of the rays going from ±ooe 17r / 4 to and from to 
±ooe-^/ 4 . 

In the setting of Theorem 3.3 we take X = R, fx the Lebesgue measure, and for ti < tj 
define 

Wu, tj = e^~^ A , K u (x,y) = K$Jx,y) := K^ c {t u x-t u y), 
W tj , ti K t .(x, y) = K^ c (tj,x; t h y). 



MULTIPLICATIVE FUNCTIONALS ON ENSEMBLES OF NON-INTERSECTING PATHS 24 

The semigroup property is obviously satisfied, while for i < j 



4 a 

OO ^ (z-z^ -|_ /• /-lOO g--V+*J 2 V° 1 



oo 



(2vri) 5 



1 /• T 00 e~ 2 r +t i ! £~ yv 1 t , 

/ / dv — - 2 = K^ c (t u x; t j ,y)+W tl ,t ] = K i^ta 3 v), 

JC J-ioo e - — +U^--xuV — U 

where the second equality follows from computing a simple Gaussian integral and the 
last equality is obtained similarly. Likewise one can check that for i < j we have 
yVti,tjYVtj,tiKp TC = Kp TC - Hence Assumption 2 is satisfied, and from Theorem 3.3 we 
deduce the following: 

Corollary 4.7. For any t± < ti < ■ ■ ■ < t n and functions qt i : R —> R, 1 < i < n so that 
Assumptions 1 and 3 are satisfied, we have 



E 



^{i-K^+Q^-^Qt, ■ ■ ■ el^-*»-i) A Q tn el^-^ A ^ rc ) i - 2( 

(44) 



In particular, the formula holds for the case qt i (x) 



L x<a, 



We do not attempt here to provide more general conditions on the functions qt t so that 
the formula holds. 

4.4. The Airyi and Airy2^i processes. All of the examples considered thus far have 
involved probability measures on ensembles of non-intersecting paths or their scaling limits. 
Going back to the discrete setting of Theorem 2.2, there was no condition that the measure 
on non-intersecting paths be positive. This condition is not met, for example, in the case 
of the Airyi arid Airy2_>i processes. These are real valued stochastic processes which are 
the scaling limits of marginals of measures (not entirely positive) on non-intersecting paths 
[9-11,42]. Even though the ensemble measure is not entirely positive, the marginal is a 
probability measure. 

We will focus on the Airy2_>.i process obtained in [11], since a similar result to that 
which we now state has already shown up in [40]. The Airy2_>.i process is a continuous 
time (non-stationary) real valued process Airy 2 _ s>1 : R — > R given by its finite-dimensional 
distributions 

p (n {Airy ^ i( ' fe) - xk} ) = det ( /_x ^ i Wi,..,MxK) ( 4s ) 

for *!<•••< t n , where xf(U, %) = ^x>xj(x) and 
(s, x; t,y) = 1 : exp ( 



VM* - s) V 4 (* 



H>s 



w 3 /3+tw 2 —yw 2yj 



l / I e 

with 

X = X-(s~) 2 , y = y-(t-) 2 , 

notation r~ = min{0,r}, and the paths 7+,7- satisfying —7+ Q j_ with 
7+ : e^+oo -)• e'^+oo, 7_ : e~'^-oo ->■ for some (f> + £ (tt/3,tt/2), ^_ £ (vr/2,7r - 

(j>+). The Airy2_s.i process crosses over between the Airy2 and the Airyi processes in the 
sense that Airy 2 _ !> i(t + r) converges to 2 1 / 3 Airyi (2~ 2 / 3 i) as r — > 00 and to Airy 2 (l;t) (the 
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Airy2 process, i.e. the top line of the multi-layer Airy2 process) when r — > — oo (in the 
sense of finite dimensional distributions). It is expected to govern the asymptotic spatial 
fluctuations in random growth models when the initial conditions are deterministic near 
the point where the hydrodynamic profile changes from flat to curved. In particular, it is 
shown in [11] that it governs the asymptotic fluctuations near the profile switch point for 
the totally asymmetric simple exclusion process starting with particles only at the even 
negative integers. 

We take again X = R and ji the Lebesgue measure, and for i < j we define 



W u , tj (x,y) = e^~ u)A (x - (t~) 2 ,y- (tj) 2 ), K ti (x,y) = K^{x,y) := K^fax^y) 

W tj , ti K ti (x,y) = Kl* 1 (t j ,x;t i ,y). 



Proceeding as in Section 4.3 one checks that these choices satisfy Assumption 2, and hence 
(under the additional assumptions) we may apply Theorem 3.3. Using the translation 
invariance of the heat kernel to rearrange the shifts appearing in the resulting formula we 
get: 

Corollary 4.8. For any t± < t-2, < ■ ■ ■ < t n , we have 



where Xi = xt - (t i ) 2 , P a f(x) = l x < a f(x) and K^^x^y) = K^ 1 (t 1 ,x + (i x ) 2 ;h,y + 



In the formula, e^'^^K^Ui should be interpreted as K|^i(*n, x + (t~) 2 ; h, y + (t x ) 2 ). 
One can use this formula directly to recover the analogous path-integral kernel formulas 
for the Airyi and Airy2 processes in the appropriate limits, and thus show that Airy 2 _ > i 
interpolates between these two processes. 

4.5. Markov processes on partitions and z-measures. The z-measures are a remark- 
able family of probability distributions on partitions that arise in representation theory 
of the infinite-symmetric group. They can be viewed as determinantal point processes on 
the one-dimensional lattice with infinite many particles, and they degenerate to a variety 
of well-known discrete and continuous determinantal point processes, see [12,13,35] and 
references therein. 

In [15], a Markov process on partitions that preserves the z-measures was constructed. 
Its dynamical correlation functions are determinantal, and they can be described via the 
corresponding extended kernel, see Section 6 of [15]. One particular limit of this Markov 
process can be seen as 'space-like' space-time sections of the multilayer polynuclear growth 
process of [36], see [14]. 

Note that it is not known how to obtain the z-measures and the corresponding Markov 
processes as a limit of an ensemble of nonintersecting paths. However, these objects can 
be viewed as an anaytic continuation of an ensemble of nonintersecting birth-and-death 
processes in the number of paths, see Section 6.5 of [15]. 

By encoding a partition A = (Ai > A2 > • • • ) by the point configuration {\ — i + |}i>i, 



the Markov process can be written as {Z(j;t) : j G Z>i,t E R}. Here Z(j;t) takes values 
in U = Z + \. 




det(J - i^i + ^ie (<2_ * l)A Ps 2 • • • e^-^-^P^e^-^K^) 



L 2 (M) 



(47) 



(*r) 2 )- 
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Similarly to (34), and with the functional q given in Definition 4.1 and operator Q given 
in (19), 



E 



n^c;;-)) 



det(l-QK^^) L2({ti _ 



t n }xZ') 



(48) 



where Kf^, g is the extended hypergeometric kernel which we will now define. 

For parameters z, z' G C such that either z' = z£C\Zorm< < m+1 for a 
m G Z, and £ G (0, 1) define a second order difference operator D zz i^ on U , depending on 
(z,z',£) and acting on functions /(•) G £ 2 (Z') as follows 



(D ZjZ ,^f)(x) = j£(z + x + \){z' + x + \)f{x + 1) 



+ yj£(z + x-\){z> + x- \)f{x - 1) - (x + £(z + z' + x))f(x). (49) 
This is a self-adjoint operator with discrete simple spectrum (1 — £)Z'. Its eigenfunctions 

tpa, 

D z , z >£lp a = (1 - £)aip a , 

are explicitly written through the Gauss hypergeometric function (see [15], equation (5.1)). 
We normalize them by the condition \\ipa \\e2(z>) = 1- Then 



K e z ^^{s,x-t,y) = < 



where Z' ± = {±± ±|, ±§, . . .}. 



(50) 



Let 0^ 



2> -"-2' 2' 

1 — ^)~ 1 D 2 2 ' j ^. Then in the setting of Theorem 3.3, (50) corresponds to 
taking X = Z', fj, the counting measure and, for i < j, 

K ti (x,y) = K z ^(x,y) := K^^0,x;0,y) = £ Mx)My), 



W tj , ti K ti (x,y) = K^^t^x-t u y) = £ er a ^- u H a {x)^ a {y). 

aeZ' + 

Thus from Theorem 3.3 we deduce the following: 

Corollary 4.9. For any t\ < t-2 < • • • < t n and functions : U — > R, 1 < i < n, so that 
Assumptions 1 and 3 are satisfied, we have 



E 



i=i 

det(J - K ZiZ , A + Q tl e (t2 - tl)D '^Q t2 ■ ■ ■ e^-^U^J^^U 



(51) 



Let us remark that if all the functions q^ have finite support then all the needed analytic 
assumptions are automatically satisfied because we are working in an L 2 space on a finite 
set. Of course, such a restriction is unnecessarily harsh, but we will not pursue this issue 
here any further. 
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Appendix A. Proofs of the results from Section 4 

We recall the following facts about trace class and Hilbert-Schmidt norms (see e.g. [43]) 
of operators in L 2 (X) for some measurable space (X, £,//), which we will use repeatedly 
without reference: 



||AB||i < p|| 2 ||5|| 2 , 
and if A has integral kernel A(x,y), 



\AB\\! < P|| op ||B||i, 



\A\ 



H{dx)n{dy) \A(x,y)\- 



1/2 



for each A,B in the appropriate space, where || • || op denotes the operator norm in L 2 (X). 

Throughout this section c and d will denote positive constants whose value may change 
from line to line. 

Proof of Corollary 4-2. Checking Assumption 2 is straightforward. We will take in this 
case Vt t = V(. = Ut i = U[. = I, and thus Assumption 1 is contained in Assumption 3, 
which we check next. 

Condition (i) is trivial. Given functions ip\ and ip2 write ipi®ijj2 for the kernel ^pi(x)ip2(y) 
and let <f> be any function with fcft 2 = 1. Then we can write QtiPk^Pk) = {Qt-fk® 
ifk), so that WQt iPk ® Wfe)||i — \\QtVk ® ^Ibll^ <8> Vfclb (note that we need to consider the 
operators with the bars because of the remark following the statement of the corollary). 
Now slightly abusing notation to write || • || 2 both for the Hilbert-Schmidt norm of operators 
in L 2 (R) and for the norm of this last space, we have 

WQufk^h = \\0--qt i )<pkh\\<l>\\2 = \\0--qu)<Pkh and lly&^lh = Ikfclhll^lh = i, 

so for 1 < i, j < n we have 



N-l 



n-1 



fe=0 k=0 

since |y>fc(x)| < cx k e~ x2 / 2 and \qti(x)\ < ce Kx2 where k < tanh(r/\/2) < \- Hence the 
only thing left to check in (ii) is that ||(5t.e _ ^ J ' _i< ^||i < oo for i < j. To that end we use 
the Feynman-Kac representation to write (setting t = h(tj — 



b(0)=x,b{t)=y 



(52) 



where b(s) denotes a standard Brownian motion and the subscript in the expectation means 
that it is conditioned (in the sense of a Brownian bridge) to go from x at time to y at 
time t. Then 

e-s/cfato 3 - 1 )*! 2 



17=7 „-£D||2 
\2 



\\Qu e ~ 



dxdyqt.ixf — e-^-y^^E 
j/w,v ) 27rt 



b(o)=x,b(t)=y 



Jinx' 



< I dxdy ce 

„2KX 2 +t 



_J--(x-y) 2 I2t w 

—e E m=xm=y 



-f*(b(s) 2 -l)ds 



< 



dx c — -=■ 



E, 



b(0)=x 



f*b(sfds 



<c' 



dx e 



2KX 2 -tanh(v / 2t)x 2 /v / 2 



< OO 



by our assumption on k, where in the last inequality we have used (1.9.3) in [18]. In the 
same way we have ||e~* D ||2 < oo and then ||Q ti e _ ( t; '~*^' D ||i < ||Qt i e~ tD ||2||e~ tD ||2 < oo. 

Finally, Assumption 3 (hi) follows from Assumption 3 (ii) thanks to the observation in 
Remark 3.2. □ 
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Proof of Proposition 4-3. Using the notation introduced before the statement of the result, 
it is clear that the functions satisfy the assumptions appearing in Corollary 4.2, so that 



E 



N n 

nii(i-%>f (**))) 
j=ii=i 



det(/-i^GUE,Ar+<5tie 



(ti-* 2 )£> . . . n P (tn-ti)D 



The left-hand side equals 



- Z1,JJ K GVE , N ) L2{]S 
(53) 



E 



N 



IJexplX) [Ml -^(ti.^f (*<))) 

3=1 



-i=l 



E 



AT 



J]exp( -5^/i(^,Af(t 4 )) + nO(<5 2 
3=1 



i=l 



-> E 



J]exp( 

3=1 



h(t,\f(t))dt 



(54) 



by the dominated convergence theorem. 

For the right-hand side of (53), writing Y h ^ = Qt^ 1 ' 1 ^ Qt 2 • • • e (t "- 1_ * n)D Q tn one 
can use the Feynman-Kac representation on each interval [*i,ti+i] as in (52) (see also (3)) 



to deduce that T 1 }'™ has kernel 

i ~-(x-y)yi(r-t) Eb{e)=xMr) __ 



where b(s) is a Brownian bridge (with diffusion coefficient 2) run from x at time t to y at 
time r. Then using (36) we deduce that 



-ph,n -ph 
1 Lr ~ 1 «,» 



1 /;(2h(s,fe( S ))+6(.) 2 -l) C 



e 2 J£ 

E? =0 Iog(l-«»(ti.6(*<)))+// 



for small enough 5. Since /i is bounded and continuous, the random variable inside the 
expectation goes to almost surely as n — > oo using a similar argument as in (54), and 
thus since this random variable is bounded by ce~ h b ^ ds the whole expected value goes 
to as n — > oo by the dominated convergence theorem. If we now define the multiplication 
operator Mf(x) = 4>{x)f{x) with 4>{x) = (1 + x 2 ) -1 / 2 then the above argument gives 
(rf,'" - T t,r) M ( x iV) ^ as n ^ oo for all x,y. To deduce that \\{T h £ ^ - V\ r )M\\ 2 -»• 
as n — > oo we use the dominated convergence theorem again together with the fact that 



t' r n - T h lr )M satisfies 



/ dxdy[(T h ^-T^ r )M(x,y)] 2 <c [ 

JM 2 JR< 



dxdye r-t E b{e)=X)b{r)=y 



e -lf[b(sfds 



< c 



< c 



< c 



dy<f>(y) 4 



1/2 



d.r 



dy e r ~ e E 



d.r 



-oo 
oo 



dye ^-t)^ m=xMr)=y 



b{l)=x,b{r)=y 
2 JT b(s) 2 ds 



-2 / / b(s) 2 ds 



1/2 



1/2 



dxE 



6(£)=x 



-2f[b(s)2 



1/2 



e - tanh (2(r-<0)x 2 



1 1/2 



< OO, 



where we have used the Cauchy-Schwartz inequality and the last equality follows from 
(1.9.3) of [18]. Checking that \\M~ 1 e^ r ~^ D Kqjj-e^ \\% < oo is simple as in the proof of 
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Corollary 4.2, so from the above we deduce that 
IK-K'gue.tv — r^V r_ ^ D ifGUE,7v) — {Kgue,n — ^£ >r e^ r ~^ D Kgue,n) II 1 

< ||(r^-r^)Af|| 2 ||M-V^K GUE ^|| 2 .0. 

ii \ i / 1 1 z 1 1 n— >oo 

Since the mapping A \-t det(I + A) L 2^ is continuous in the space of trace class operators 
(see [43]), we deduce that the right-hand side of (53) converges to det(I — -Kgue,at + 
r er e( - r ~ e)DK GUE,N)L 2 {R), and hence 



E 



N 

JJexp 



h(t,Xf(t))dt 



det(J - K G UE,iv + T h er e^-^ D K GlJE 



□ 



Proof of Corollary 4-6. Fix / G L 2 (M) and write (f){x) = e rx l x >o + tp(x) 1 l x <o- Then for 
i < j (note that as in the proof of Corollary 4.2 we need to consider the operators with 
bars), writing /(A) = dx Ai(x + X)f(x) we have 

-i 2 



\\Q t / u -^ H f Wi < 



dx 



< c 



dx 



oo 
oo 



dy dX (p(x)e 



A\{x + X)M{y + X)f{y) 



dXMxYe 



2 p 2X{t j -t i ) 



Ai(x + Xf 



dXf(X)' 



(55) 



dx^ixfe- 2 ^- 1 ^ \ dXe^ x+1 ^-^ Ai(A) 2 < c' 



where in the last line we have used the Parseval identity for the Airy transform f f 2 = f f 2 - 
The fact that d < oo follows from the assumption on r and ip and the bounds 

|Ai(s)| < ce _ i x3/2 forx>0 and |Ai(x)| < c|x| _1/4 for x < (56) 

(see (10.4.59-60) in [1]). This shows that for i < j, Qt i ^~ lytj ~ ti ^ H is a bounded operator 
mapping L 2 (IR) to itself. Similar computations allow to check the rest of Assumption 
l(i). To check (ii) we use the formula given in Lemma 3.1. Each term can be written 
as a product of the form (-l) fc (e (<1 ~ <a « )H Q tao ) • • • (e {ta ''-i~ ta ^ H Q ta J(e^ a k- t ^ H K 2 ) with 

1 < a® < ■ ■ ■ < afc < n. The k + 1 factors coming after (— l) k can be checked to be bounded 
operators on L 2 (M) by a computation similar to (55), and a simpler computation gives the 
same for e ( - ta k- t ^ H K 2 using the spectral formula for its kernel. 

As in the previous example, checking Assumption 2 is straightforward. For Assumption 
3 we choose 

VtJ(x) = U' t J\x) = ^(x)f(x), Vlj{x) = U t J(x) = ^{x)- l f{x) 

with ip(x) = e~ rx / 2 l x >o + </?(z) 1//2 l x <o- Condition (i) is obvious. Now note that K 2 = 
B P Q B , where B Q (x, A) = Ai(x + A) and P a f(x) = f(x)l x > a , so 

WVt^KzVl^ < ||V^i?oPo||2||PoA)V£||2. (57) 

We have 



V u Q t .B Q Po\\i < c dx dX^xY^xY ki{x + XY 



f'OO f'OO f'O f'OO 

■ dx dXe rx Ai(x + A) 2 + c dxifix)^ 1 dX Ai(A) 2 . 

JO Jo -'-oo J -x 



The first integral on the right-hand side is clearly finite by (56). For the second one, note 



that by (56) f~ dX Ai(A) 2 < c(l + 



|l/2> 



so the integral is finite by the assumption on ip, 
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which by the Cauchy-Schwartz inequality implies f° dx \x\ 1 ^ 2 Lp{x)~ l < oo. H-Ffci-BoVt- Ih < 
00 follows from the exact same calculation, and hence from (57) we get that V^Q^-K^V^. 
is trace class. The same proof shows that Vt i Qf .e^'~ ti J H K2VI. is trace class for i < j. To 
check that Vt i Q t .e~ tH V^. is trace class for i < j and t = tj — t$ > we start by writing 

Wt^e-^h < \\V tl Q u e- tH /%\\e-^V>\\ 2 . 
For the first factor we use again the explicit formula for the kernel of e~ tH to obtain 

\\V ti Q u e~ tH/2 \\2 = [ dx dy dXdaij(x) 2 <P{x) 2 e t( - x+a ^ 2 Ai{x+X)A[{y+X)Ai(x+a)Ai(y+a) 

dxiP{xf<t>{x) 2 e- tx / dX e <A Ai(A) 2 

-00 J —00 



00 

dx 

-00 



(58) 

which is finite by the similar arguments as above. \\e~ tH l 2 V(. H2 can be bounded in the 
same manner, and we deduce that V^Qt e~ tH V^ is trace class. The same proof works for 
V ti Q ti e- tH K 2 V{.. 

To get (iii) we use Lemma 3.1 and rewrite each term in the sum as 

(-l) fc (L/ ti e^- i -)^ ao ^J(^ oe (^o-^ 1 )^ ai ^ i) . . . (U^/^-^Q^UiJ 

Since Ut t = V{. and Up = Vj^, each factor above corresponds to the adjoint of one of the 
factors appearing in (ii). Since the adjoint of a trace class operator is also trace class, we 
deduce that the whole product is trace class. □ 

Proof of Corollary We already indicated how to check Assumption 2. One checks 
directly that the first three operators in Assumption l(i) are bounded operators preserving 
L 2 (R), while the last one can be checked using (46) and arguing about the Airy functions 
appearing there similarly as in the previous proof. Assumption 1 (ii) follows similarly using 
Lemma 3.1. Assumption 3 can be checked following the same ideas as in the proof of 
Corollary 4.6 and using the arguments in Appendix A of [9] to provide the necessary 
analytical estimates. □ 
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